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Abstract 



Band Husimi distributions (BHDs) are introduced in the quantum-chaos 
problem on a toral phase space. In the framework of this phase space, a quan- 
tum state must satisfy Bloch boundary conditions (BCs) on a torus and the 
spectrum consists of a finite number of levels for given BCs. As the BCs are 
varied, a level broadens into a band. The BHD for a band is defined as the 
uniform average of the Husimi distributions for all the eigenstates in the band. 
The generalized BHD for a set of adjacent bands is the average of the BHDs 
associated with these bands. BHDs are shown to be closer, in several aspects, 
to classical distributions than Husimi distributions for individual eigenstates. 
The generalized BHD for two adjacent bands is shown to be approximately 
conserved in the passage through a degeneracy between the bands as a non- 
integr ability parameter is varied. Finally, it is shown how generalized BHDs 
can be defined so as to achieve physical continuity under small variations of 
the scaled Planck constant. A generalization of the topological (Chern-index) 
characterization of the classical-quantum correspondence is then obtained. 

PACS numbers: 05.45.+b, 03.65.Sq 
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I. INTRODUCTION 

The main objective of "quantum chaos" is to understand the correspondence between 
classically nonintegrable systems and their quantum counterparts in the semiclassical limit 

During the last two decades, significant progress has been made in the study of this corre- 
spondence with the discovery of phenomena such as dynamical localization "scarring" 
of eigenstates by unstable periodic orbits , and statistical properties of the eigenspec- 

trum [T|,0] . However, the relation between classical phase-space structures and corresponding 
quantum-dynamical entities is still far from being completely understood. 



In this paper, we introduce quantum-mechanical distributions which, in the semiclassical 
limit, are expected to approach in a natural way classical distributions on both regular and 
chaotic phase-space structures. This will be done in the framework of a toral phase space, 
where, as shown in a recent series of works [p|-[TO[|, some interesting new insights in the 
quantum-chaos problem can be achieved. Quantum dynamics can be reduced to a torus if 
two conditions are satisfied (see Sec. II for more details). First, the classical map for the 
system is strictly periodic in all the phase-space coordinates. The simplest nonintegrable 
system possessing this property is the "kicked-Harper" (KH) model |[^-p!5| with Hamiltonian: 

oo 

H = Acos{v) + Acos{u) ^ 5{t/T - s) , (1) 

S = — OD 

where u and v are dimensionless conjugate phase-space variables (with Poisson bracket 
{u, v} = 1/1,1 being some classical action), A is the amplitude, and r is the time period. 
The system (P is exactly related to the problem of periodically kicked charges in a 



uniform magnetic field under resonance conditions [|T6[. In the limit r 0, (|l|) reduces to 
the integrable Harper Hamiltonian Hq = A cos{u) + A cos{v) . The transition from integrable 
to chaotic phase-space structure, as the dimensionless classical parameter 7 = At/{27tI) is 
increased from 0, is shown in Fig. 1. The second condition for quantum dynamics on a 
torus is that a scaled h, denoted here by p {[u, v] = 2Tcip = ih/I), assumes rational values, 
p = q/p {q and p are coprime integers). The admissible quantum states are then those which 



satisfy Bloch quasiperiodic boundary conditions (BCs) on the torus, see Sec. II. The energy 
or quasienergy spectrum consists of precisely p levels, and, as the BCs are varied, each of 
these levels spans a band. 



The advantage of this framework is that it allows for a characterization of the classical 
quantum correspondence by means of integer topological invariants, the Chern indices UsHlO 



associated with the p bands. The Chern index a for a band is analogous to the quantum 
Hall conductance carried by a magnetic band in a perfect crystal ||T7|-p4|, and is a measure 



of the sensitivity of the eigenstates in the band to variations in the BCs [pl- [l0|j25|] . For g = 1, 
the toral phase space coincides with the basic unit cell of periodicity of the system. In this 
case, where the classical-quantum correspondence can be established in the simplest and 
most natural way, a can assume, in principle, all values. Several arguments [3,p|,p5| , sup- 
ported by numerical evidence, then indicate that if the Husimi distribution of an eigenstate 
is localized, in a semiclassical regime (p <^ 1), on classical regular orbits [e.g., Kol'mogorov- 
Arnol'd-Moser (KAM) tori or periodic orbits] the corresponding band has a = 0. On the 
other hand, eigenstates whose Husimi distribution is spread over the classical chaotic region 
should correspond to bands with cr 7^ 0. The transition from a nearly-integrable regime [e.g.. 
Fig. 1(a)], where almost all a = 0, to a fully chaotic regime [e.g.. Fig. 1(d)], where almost 
all cr 7^ 0, as a nonintegrability parameter is varied, takes place via degeneracies between 
adjacent bands. In the passage through a degeneracy between bands h and 6', the Chern 
indices cif, and ay change, respectively, by ±Acr, where, generically [|1^,^], |Acr| = 1. A 
"diffusion" of Chern indices [p^ occurs then in the transition above. 



Despite this characterization of the classical-quantum correspondence for small p, the 
eigenstates may not be considered close to classical phase-space structures, strictly speak- 
ing, for any p. This is because of the following reasons: (a) While the BCs satisfied by 
a quantum state have a well-defined physical meaning (see Sec. HI), they are of a purely 
quantum nature. In particular, the strong dependence of the eigenstates on the BCs for 
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cr 7^ 0, and, in several cases (see Sec. Ill), also for a = 0, has no classical counterpart, (b) 
The Husimi distribution of an eigenstate always assumes p zeros in the torus, see Sec. II. Be- 
cause of this fact, an eigenstate cannot tend, in the semiclassical limit, to the microcanonical 
uniform distribution on the chaotic region P,p8|. (c) In the general case of g 7^ 1, the exact 
eigenstates may be viewed as arising from quantum tunnelling between degenerate classical 
orbits located in q adjacent unit cells iPHTlI]. As a consequence, a is always nonvanishing 
T0| and the topological characterization of the classical-quantum correspondence cannot be 



extended straightforwardly to this general case. 



We show in this paper that a natural way to overcome these difficulties is simply to 
average uniformly the Husimi distributions of the eigenstates in a band over all the BCs. In 
what follows, we refer to the result of this averaging as the band Husimi distribution (BHD). 
The BHDs turn out to be closer, in several aspects, to classical distributions than Husimi 
distributions for individual eigenstates. It is well known that smoothing a probability 
distribution over a range AE of energy levels has the effect of washing out purely quantum 
structures such as scars. As AE is increased, this effect increases and the smoothed proba- 
bility distribution (e.g., the spectral Wigner function in Ref. [Q) becomes closer to a classical 
distribution. In our case the smoothing is done not over a range of discrete energy levels 
but over the continuous range of one band, corresponding essentially to a single level in the 
framework of a toral phase space. This is thus the minimal smoothing in this framework, 
and it is performed just for the sake of eliminating the purely quantum effects of individual 
BCs. The BHD may be viewed as the representative Husimi distribution for a level in the 
torus. In several important cases, we shall find it necessary to generalize the concept of 
BHD by considering the average of BHDs associated with a set of adjacent bands. 



This paper is organized as follows. In Sec. II, we summarize the relevant known facts 



about quantum dynamics on a torus PHT3- Sec. HI, the concept of BHD for a single 



band is introduced. In Sec. IV, the concepts of generalized band (set of adjacent bands) 



and the associated BHD are introduced and studied. In Sec. V, the generahzed BHD for 
two adjacent bands near a degeneracy of these bands is studied. In Sec. VI, we show how 
to define generalized BHDs for g 7^ 1 in order to achieve physical continuity under small 
variations in p. In this way, the topological characterization of the classical-quantum corre- 
spondence is extended to the general case of g 7^ 1. Conclusions are presented in Sec. VII. 



II. QUANTUM DYNAMICS ON A TORUS 



We provide here some relevant background on quantum dynamics on the torus [§-[T0 
using, for later convenience, the notation of Ref. ||10|, where a general formulation of the 



problem was presented. Consider a classical area-preserving map strictly periodic in the 
phase space (u, v) with a 27rx27r unit cell, which is the basic torus T^. The one-step quantum 
evolution operator corresponding to the classical map is U{u, v), where [u, v] = 27rzp. If p 
is a rational number, p = q/p {q and p are coprime integers), there exists a pair of "smallest" 
commuting phase-space translations 

Di = e^lf, 1)2 = e'P^ . (2) 

Because of [u, v] = 27rip, Di (-D2) is a translation by 27r (27rg) in the v (u) direction. Since 
U{u, v) is periodic in both u and v with period 2tt, it commutes with Di and D2- There 
exist therefore simultaneous eigenstates \1' of U, Di, and D2- 



^i|^6,w) = e™i/''|^,,w) , ^2|^6,w) = e'P^'l^b,^) , (3) 

where 6 is a "band" index (see below), u;f,(w) is the quasienergy, and w = {wi, W2) is a 
Bloch wavevector varying in the "Brillouin zone" (BZ) < Wi < 27rp, < W2 < ^n/p. 
Now, for each given value of w, Eqs. @ can be interpreted as quasiperiodic boundary con- 
ditions (BCs) satisfied by the eigenstates \'^b,v/) in the "quantum" toral phase space Tq. 
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< M < 27rg, < V < 27t. It can be shown that, for each w, the quasienergy spectrum 
consists precisely of p levels Ubi^), b = 1, p. As the BCs are varied (by varying w in the 
BZ), each level broadens into a "band". 



In the absence of band degeneracies [oobiyv) 7^ t^6'(w) for all w and 6' 7^ 6], \'^b,w) must 
be periodic in the BZ up to a constant phase factor depending, in general, on w. The phase 
of l^^ft.w) may be chosen so that \'^b,w) will be exactly periodic in one direction, say wi, but 
then it will be periodic in W2 only up to a phase factor which, in its simplest form, has a 
phase that is linear in wi ||10|| : 

\'^b,wi+277p,W2) = \^b,w) , (4) 

|*6,«,l,«'2+27r/p) = exp(i(Tfew;i/p)|^b,w) • (5) 

Here the constant ab must be an integer in order for Eq. to be consistent with Eq. 
(|). It is easy to see from Eqs. (|^) and (|) that 2TTab is the total phase change of \'^b,w) 
when going around the BZ boundary counterclockwise. This phase change is independent on 
phase transformations of \'^b,w)- In fact, the integer ab is a topological number, the Chern 
index, which can be expressed in a form manifestly invariant under phase transformations 
[see Ref. and expression (^) below]. 



The general form of |^E'b,w) m the ^-representation is [|T0 | 



(6) 



m=0 



where 0b(m; w) are expansion coefficients and ijj-wiv) are A;g-functions 



^w('y) = E ew{i^uii/q)6{v - W2 + 27rl/p) . 

l=—ao 



(7) 



Without the label b, (^) is the most general expression of a quantum state in the v- 
representation for given BCs (i.e., given w). A completely analogous expression is obtained 
in the conjugate {u) representation by Fourier-transforming (P). These expressions imply, 



essentially, a discretization of the toral phase space Tq into points, Umi = wi + 27rmip, 
= + 2711712/ p, ^1, m2 = 0, p — 1, with 4>{mi; w) being the probability amplitude 
for u to assume the value Um^. The corresponding amplitude for v = is the discrete 
Fourier transform of (f){mi] w). This discretization was first studied by Hannay and Berry 
Pm in the special case of g = 1 (Tq = T^) and w = (strict periodicity). They used the p 



amplitudes 0(m; w = 0) to investigate the evolution of a quantum state in the cat maps. 

A most useful representation of |^'6,w) is the coherent-state representation \E'fe_w(R') = 
(R|^fe_w), where R = (m, v) and |Ro) is a coherent state. 



2 \ 1/4 



[v - vo/2) 



(8) 



47rp 27rp 

Here a is the "squeezing parameter" , related to the parameters of the harmonic oscillator for 
which the quantum evolution of the state (H) is nondispersive. An equivalent expression for 
a coherent state depends on uq and Vq only through the complex number zq = uo/a — iavo 
(see, e.g., Refs. [^,^), giving an analytic representation "^b^^i^z). For future convenience, 
however, we shall use the representation \E'{, w(R)7 choosing, as in Ref. PB|, the symmetric 
value a = 1. 

Important properties of \&f,^w(R) concern its zeros, in both the R and w variables. At 
fixed w, w(R) always assumes exactly p zeros R = Roj(w) (j = 1, p) in Tq (counting 
possible but nongeneric zero multiplicities). This property was proven in Ref. for q = 1 



and a = 1, but it can be easily generalized to all q and a, using Eq. (|]), together with 
Eqs. d^) and (§). It is also well known that the p zeros Ro,j(w), like the p amplitudes 
(pbiiTT-', w) in (0), completely determine the wavefunction ^{, w(R)- At fixed R, the number 



iVo(R) of zeros w = wo(R) of \l';,^w(R) in the BZ is not smaller than |crf,| |3^. These prop- 
erties imply that the p zeros Roj(w) must wind around Tq at least |crb| times when w is 
varied in the BZ. Thus, if for some R = R' \l'5^w(R') does not vanish in the BZ, ab = 0. We 
shall refer to the union of all R' as the localization domain of \E'5^w(R), for reasons that will 
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become clear in the next paragraph. As w is varied over the entire BZ, the p zeros Roj(w) 
never enter this domain (see also Ref. p5|). 



In a nearly- integrable situation [e.g., very small A in the KH model (|1])] and in a semiclas- 
sical regime (p ^ 1), most eigenstates ^^b^wlR-) will be localized on regular classical orbits 
and the corresponding bands are usually very narrow (almost independent of w). A good 
approximation to the Husimi probability distribution |\E'b^w(R')P in its localization region is 



given by ||8|,|34 



|*6,w(R)r ~ T^^exp ^ 



(9) 



V{R) 1 2np 

where Ho(R) is the energy function for an effective integrable Hamiltonian Ho in this nearly- 
integrable situation, Ei, is an energy eigenvalue of Hq approximating the very narrow band 
[Eb ^ uJb{'w)h/T], V(Il) is the phase-space velocity for Hq, and A/" is a normalization con- 
stant. Rel. d^) manifestly shows that the eigenstates are quite insensitive to variations in 
w in the region of phase space where they are localized. In particular, the p zeros Ro,j(w) 
should never enter this region, thus implying a finite localization domain (as defined above). 
We therefore expect that eigenstates that are localized on regular classical orbits for p ^ 1 
should belong to bands with ai, = (see also Sec. III). 

It is important to remark here that the value = may occur only for q = 1, since ab 
has to satisfy the Diophantine equation pab + qfib = 1, where pb is a second integer |TO|,pn 



For q ^ 1, the exact eigenstates are quite "nonclassical" , since they may be viewed as arising 
from quantum tunnelling between q degenerate classical orbits located in the q adjacent unit 
cells defining Tq p|-|TT|. Thus, even in a completely integrable situation, the p zeros Roj(w) 



always cover the torus Tq when w varies in the BZ 0. 

In a strongly chaotic situation, where a typical orbit fills the entire phase space T^, the p 
zeros Roj(w) for q = 1 are expected to be distributed almost uniformly in Tq = and 



to explore the entire as w is varied in the BZ pj. Almost all the bands should therefore be 
characterized by nonvanishing at- The transition from a nearly- integrable situation (almost 
all o"6 = 0) to a strongly chaotic one (almost all o";, 7^ 0), as a nonintegrability parameter 
is varied, takes place via successive degeneracies between adjacent bands [§,0 (see also the 
Introduction). 



III. BAND HUSIMI DISTRIBUTIONS 

We start this section by considering in more detail the nature of the BCs (§) using the 
KH system (|I]) as a model. The nearly-integrable regime for this system corresponds to very 
small values of the classical parameter 7 = At/{27iI) [see Fig. 1(a)]. In this regime, the 
system is well described by the Harper Hamiltonian Hq = Acos{u) +Acos{v). For p = 1/p, 
p odd, the energy spectrum of Hq consists of p bands Eb{w), 6 = 1,..., p, in order of increas- 
ing energy. Only the central band [b = {p + l)/2] has a nonvanishing Chern index cxb = 1 
P^ . In the semiclassical regime, p ^ 1, the eigenstates of this band are concentrated on 
the separatrix orbit [see Fig. 1(a)], which is not contractible to a point. On the other hand, 
the eigenstates of the other bands, that have vanishing Chern indices, are concentrated on 
orbits which are contractible to a point. 



This state of matters persists also for 7 not very small, when the separatrix orbit breaks 
into a stochastic layer, and island chains emerge from the Hq contractible orbits. For exam- 
ple, for p = 1/11, the Chern indices appear to be the same as in the Harper case in the 
entire interval 7 < 0.26 P,|35[[! Namely, only the central quasienergy band features a nonva- 



nishing Chern index = 1- We have studied numerically the sensitivity of the eigenstates 
to variations in the BCs (i.e., variations in w) for several values of 7 in the interval above. 
We find that eigenstates in o";, = bands sufficiently "far" from the central band (i.e., b 
close to 1 or 11) are indeed almost insensitive to variations in w (large localization domain). 
This is not the case, however, for bands sufficiently close to the central band. While the 

10 



Chern index for these bands vanishes, the sensitivity of the eigenstates to variations in w 
is quite strong (very small or empty localization domain), almost as that of eigenstates in 
the central band. This is clearly illustrated in Figs. 2 and 3. In general, we expect strong 
sensitivity to variations in w in o";, = bands that are sufficiently close to CTf, 7^ bands. 



This sensitivity has no classical analogue, since w is a purely quantum characterization 
of the eigenstates, which are concentrated on different regions of Tq for different values of 
w. In fact, by comparing Eqs. @) with the definitions @) of the operators IDi and 1)2, we 
see immediately that w is just a quasicoordinate of R = (u, v): 

Wi <^==^ pi— mod 27r I , W2 <^==^ -(pv mod 2n) . 
\P J P 

To understand this better, consider, for given w = w', the "displaced" quantum state 

^(2:o)|^6,w') = exp{zoa^ - z^a) \<i/b,^>) , (10) 

where Zq = Uq + Wq and a = {u + iv)/ {4:TTp) is the annihilation operator. It is easy to verify, 
using the well known commutation relation for the phase-space translations D{zo) that 
the state (p!0|) satisfies the BCs with the displaced value of w = w' + Rq [Rq = {uq, vq)]. 



On the other hand, using the expression (H) (for a = 1), we find that the coherent-state 
representation of the state ([ToD is given by 

(R|^(2;o)|^b,w') = exp[i(Mt;o - Uov)/{A7Tp)] ^6,w'(R - ^o) • (H) 

Thus, the shifted function \E'b ^'(I^ ~ I^o) iS; up to a phase factor, the coherent-state repre- 
sentation of a quantum state characterized by the displaced value of w = w' + Rq. This is 
a vivid illustration of the notion of "quasicoordinate of R" for w. 



From Eq. ( [Tl| ) one may get the impression that there is always strong sensitivity to 
variations in the BCs. However, since D{zq) does not commute, for general zo, with the 
evolution operator 0, the displaced state (|r^) will not be, in general, an eigenstate. Still, 
one has the expansion 
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p 

b'=l 

where the coefficients Cft_fe/(Ro, w) = (\&{,/ w|-D(2;o)|^b,w')- Assume, for example, that \E'b^w(R') 
is almost insensitive to variations in w (large localization domain). Then, from Eq. (pH]) and 
w = w' + Rq it follows that the coefficient Cb^b(R-o, w) should be relatively small for almost 
all w if Rq is larger than the typical width of the localization domain. If, on the other hand, 
\E'f,,w(R.) is highly sensitive to variations in w, the coefficients C6,fe/(Ro, w), for b' ~ b, may 
be relatively large for "many" pairs (Rq, w) 7^ 0. For example, in the cases shown in Figs. 
2 and 3, these pairs include [Rq = (vr, tt), w = (0, tt/II)] [i.e., w' = (vr/ll, 0); compare 
Fig. 2(b) with Fig. 2(c) and Fig. 3(b) with Fig. 3(c)] and [Rq = (tt, tt), w = (tt/II, tt/II)] 
[i.e., w' = (0, 0); compare Fig. 2(a) with Fig. 2(d)]. 



Because of the purely quantum nature of w and the BCs, it is natural to perform a 
uniform average over w in the Brillouin zone in order to obtain "more classical" quantities. 
In this paper, we shall consider the average of the Husimi probability distribution |\I'5^w(R)P, 
giving the band Husimi distribution (BHD) for band b: 

P,{R) = -l- f dw\%,^{R)\\ (12) 
|B/| Jbz 

where |BZ| = An'^q/p'^ is the area of the Brillouin zone. The BHD ([T^) corresponds to the 
minimal smoothing of a probability distribution in the framework of a toral phase space, 
namely the smoothing over the continuous range of a single band. We now show that 
this smoothing is sufficient to make the BHD "more classical" than an individual Husimi 
distribution |\I';,vv(R)P in several aspects. First, we notice from Eq. that |\E'f,_w(R)P 
is periodic only with unit cell Tq, i.e., the quantum phase space, which differs from the 
classical one, T^, whenever q ^ 1. Now, using the relation [rU 

Z)(-27rs)|^5,w) = exp{isw2/ p)\'^b,w,~27Ts,w2) 
[s integer) as well as Eq. ( [TI]) in ([12|) , we easily find that 
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P,{R) = -J2Pi'\u + 2ns, v) , (13) 

1 s=0 

where P/'^^(R) is defined as in Eq. (]T^), but the integral is performed over 1/g of the BZ, 
i.e., < wi, W2 < 27r/p, and |BZ| is replaced by |BZ|/g. It is then clear from Eq. (p!3D 
that -Pb(R), unlike |\l/fc,w(R-)P, is periodic with unit cell for general q. This allows one to 
impose on Pb(R) the normalization condition 

/ rfRPb(R) = 1 . (14) 

This makes P6(R) analogous to a classical probability distribution on the classical phase 
space T^. 



Second, |\E'5w(R)P always assumes p zeros Roj(w) (j = 1,..., p) in Tq (see Sec. II). 
These zeros give |\I'5vv(R)P a rather "nonclassical" appearance, for example, they do not 
allow |\I'b,w(R)P to approach, in the semiclassical limit, the microcanonical uniform dis- 
tribution in the chaotic region (strong-chaos regime) P,^. On the other hand, the BHD 



never vanishes in the phase space [P6(R) > for all R in T^], simply because the p zeros 
Roj(w) (j = 1,..., p) generally vary with w, and, by definition [see Eq. (|12])], a BHD 
involves an integration over all w. It is therefore possible for a BHD to resemble a classical 
probability distribution. For example, in Fig. 4 we show the BHD for the case consid- 
ered in Fig. 2 (6 = 6). The relatively high probability density near the hyperbolic points 
R = (tt, 0), (0, tt) can be easily understood from classical considerations. The band b = 6 
corresponds to a "broken" separatrix orbit (that is, a homoclinic orbit in the chaotic layer), 
and the phase-space velocity on this orbit vanishes as one approaches the hyperbolic points. 
Accordingly, the approximate formula (^ suggests that the BHD should assume relatively 
high values near these points. Fig. 5 shows that the localization region of the BHD for 
band 6 = 5 is completely different from that of individual Husimi distributions (see Fig. 
3). Although this band has 0"^ = 0, it exhibits strong sensitivity to variations of the BCs. 
The representative or dominant localization region for such a band can be found only by 

13 



inspecting its BHD. In the strongly chaotic regime (7 ^ 1) and in the semiclassical hmit, 
the BHDs are expected to approach the microcanonical uniform distribution. 



Finally, consider the special but important case of bands with at = 0, which is possible 
only for g = 1 (see Sec. II). Here, an eigenstate ^^ft^wlR-) can always be written as a 
symmetry-adapted sum [pO| , p^ 

00 

^5,w(R)= E exp[-ip{hwi + l2W2)]{R\D['&i\vb) , (15) 

h ,/2=— 00 

where Di and D2 are the basic phase-space translations (Q) and | is some square-integrable 
state, which is analogous to a Wannier function p0| , p2| , p4| . Inserting Eq. (^) into Eq. (^2|), 
we easily obtain the exact expression 

00 

Pfe(R)= \{u + 27Tli, v + 27rl2\ipb)\^ . (16) 

While the Wannier function (R| is not invariant under gauge transformations in which 
the eigenstates are multiplied by phase factors exp[i0(w)] P3|, the BHD (ITB) is gauge in- 
variant. In a nearly-integrable situation and in a semiclassical regime, |(R|v56)P may be 
identified with the "quasi-mode" of Ref. which is well localized on a classical regular 
orbit and, in the limit p ^ 0, tends point-wise to zero outside this orbit. Similarly, the 
BHD (|16D tends point-wise to zero outside the periodic repetition of the orbit on all unit 
cells (/i, 12). It is therefore a periodic version of the quasi-mode, appropriate for a toral 
phase space. A good approximation to the BHD should be given by the right-hand side of 
Eq. (P), which, like (0), is essentially independent of w and is periodic with unit cell T^. 
The difference Zb^w(R) = |^b,w(R)P — ^(R) is the sum of the overlaps of the quasi-mode 
{u + 2ttIi, V + 27Tl2\ipb) [ill unit cell (/i, ^2)] with a quasi-mode in a different unit cell. The 
function Zb w(R) is thus of a purely quantum nature and it is entirely responsible to the 
p zeros of \l/6,w(R) and to the sensitivity of \E';,,w(R) to variations of w for R outside the 
localization domain. This clarifies the classical nature of the BHD in this case. 
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IV. GENERALIZED BANDS AND BHDs 



In several important situations, some of which will be considered in the next sections, it 
is necessary to generalize the concept of BHD by smoothing over more than a single band, 
namely over a set of adjacent bands b = bi, 6jv. This gives the generalized BHD 

n,-.,(R) = ^E^^(R) • (17) 

b=bi 

The additional smoothing over bands should give a "more classical" BHD, as when smooth- 
ing over many levels in a bounded quantum system 0. The "maximal" smoothing is, of 
course, that over all the p bands. From the completeness of the eigenstates @, together 
with the normalization condition (|14|) , we find in this case that 

i|:n(R) = ^, 

where |T^| = 47r^ is the area of the classical torus. Thus, as one could expect, the generalized 
BHD in this case is just the uniform distribution over phase space. 



The set of adjacent bands can be considered as a single entity, a generalized band 
(GB). One may perform arbitrary linear combinations of GB eigenstates at fixed w to 
obtain general states satisfying given BCs on the torus. While these states are generally 
not eigenstates of the evolution operator, they are "almost stationary" provided that the 



energy or quasienergy width of the GB is sufficiently small |10|. The set of all these states, 
for all w, is the space of the GB. A natural starting basis for this space is, of course, the set 
of eigenstates \E'f,,w(R), b = bi, b^, at each fixed value of w. An arbitrary basis will 
then be given by 

^„,w(R) = E 4"V)*b,w(R) , (18) 

b=bi 



n = 1,..., N. To ensure orthonormality of the basis ( [T8| ) in the new "band index" n, 
the coefficients -B^"^(w) must form a unitary matrix. Obviously, the states satisfy the 
BCs (H). In addition, it is natural to require that these states will satisfy quasiperiodicity 
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conditions in w, analogous to those of Eqs. and (H), with well-defined Chern indices an- 
Clearly, this will be the case only if the matrix B^^^ (w) in Eq. (p!8| ) satisfies these conditions 
with Chern indices ab,n = o"„ — ab- The determinant of this matrix must be strictly periodic 
in w, otherwise it will vanish at some w (see Sec. II and note |22|), which cannot happen 
since the matrix is unitary. It follows from this that 

bjv N 
6=61 n=l 

In other words, the total Chern index ugb of the GB is "conserved" under the basis trans- 
formation in Eq. (p!8|) . 



Another "conservation" law following from Eq. (|T8]) is 

Ei^uR)r = Ei^n,w(R)r 

6=61 n=l 



(20) 



In particular, Eq. (^), when integrated over w, implies that the generalized BHD ( p!7[ ) can 



be calculated using an arbitrary basis (|TH). 



An important case is when one can find a new basis (IB) whose Chern indices ct„ all 
vanish. Before discussing the meaning of this case, we first determine the conditions that 
need to be satisfied to make it possible. Because of Eq. ([T9|) , a necessary condition is clearly 
that (Tgb = 0. This condition is also sufficient, since if (Tgb = one can always find a unitary 
matrix bI"\w) with Chern indices cTfe^„ = — cTf, (implying that = for all n). In fact, one 
can simply choose Bj^"'\'w), n = 1, A^, as the A^ orthonormal eigenvectors of a A^ x A^ 
Hermitian matrix that is strictly periodic in w and whose A^ homotopic invariants (Chern 
indices) are — at, b = bi, bj^. Such a matrix can be exphcitly constructed for any given 



set of integers cr^ with o"gb = |18]. It is worthwhile to stress here that due to the general 



Diophantine relation pat + q^ib = 1 P^ , |20| the condition ctgb = can be satisfied only if A^ 
is a multiple of q, i.e., the minimal N is N = q. 
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The fact that one can find a new basis ([ISD with ct„ = for all n means that the GB can 
be viewed as "weakly sensitive" to variations in the BCs, despite the fact that the original 
Chern indices ah may be all different from zero. This can be expressed in a more precise way 
using Eq. pO|). Since d"„ = 0, one can assume that \l/„,w(R') has a finite localization domain 
(see Sec. II). Eq. ( pOD then implies that the function Y!blbi 1^6,-^(1^)1^, characterizing the 
GB, has also a finite localization domain. In this sense, the GB exhibits "weak" sensitivity 
to variations in the BCs. Nonzero values of the Chern indices ah in this case only mask the 
true nature of the GB, which is best described in terms of the new basis. One can thus say 
that a GB consisting of q adjacent bands with (Tgb = is analogous to a band with cxfe = in 
the case of g = 1. In Sec. VI, these ideas will be further developed in order to generalize the 
Chern-index characterization of the classical-quantum correspondence to the case of g 7^ 1. 



V. BHDs NEAR DEGENERACIES 

In this section, we show that the generalized BHD for two adjacent bands (to be denoted, 
for simplicity, by 6 = 1, 2) is approximately conserved as a nonintegrability parameter 7 
is slightly varied through a degeneracy point 70 of these bands, i.e., ci;i(wo) = c<j2(wo) at 
7 = 7o, where wq is some isolated value of w. This despite the fact that the separate BHDs 
of the two bands usually change significantly under such variation. 

Formally, the unitary evolution operator can be written as f/ = exp[— ^(7(7)], where G(7) 
is a Hermitian operator. The quasienergy states (H) are eigenstates of 6(7) with eigenvalues 
uih(yf). Consider two values of 7, 71 and 72, very close to 70 and such that 7i < 7o < 72, 
and let us denote by \^^^\,), j = 1, 2, the eigenstates of 6(7^). We assume that for 7 in the 
interval [71, 72] (containing the degeneracy point), and for all w, the distance between ti;i(w) 
and ti;2(w) is significantly smaller than the distance between any of these quasienergies and 
ci;f,(w), 6 7^ 1, 2. In this case, one can write, to a good approximation. 
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*I',l(R)^E<^(w)<i(R), (21) 

6=1 



where the expansion coefficients 5^ (w), b, b' = 1, 2, form a 2 x 2 unitary matrix built from 
the normahzed eigenvectors of the 2x2 Hermitian matrix Gb,b'{^) = (^6^w|G^(72) l^i^w)- 
Eq. (plD is thus an approximate special case of Eq. (0). Nevertheless, the general relation 
(p^) holds exactly in this well. It expresses the well known conservation of the total 



Chern index ai + 02 in the passage through a degeneracy point [lT9| , p6 



From Eq. ( |21| ) we get the relation 

+ iv&w (R)p ^ \^%m^ + iM1(R)r , (22) 

which is an approximate special case of Eq. ( PP] ) . After integrating Eq. ( [2^ ) over the entire 
BZ, we obtain the approximate conservation law for the generalized BHD of the two bands: 

Pg(R)^Pg(R). 

As a ffist, instructive example, we consider the degeneracy between bands 6 = 3, 4 in 
the KH model with p = 1/11 for 7 = 7o ~ 0.264 (see Ref. |]). For 7 = 71 = 0.26 < 70, the 
Chern indices of the two bands are, respectively, cr-s = = 0, while for 7 = 72 = 0.2645 > 70 
they change to (Xs = —cr^ = 2 For 7 = 73 ^ 0.2653, the Chern indices re-assume the 



values (J3 = (74 = due to a second degeneracy between the two bands. Thus, for both 
7 = 7i and 7 = 73, one can associate Wannier functions (R\'^b) with these bands, as in Eq. 
(P^). These functions are expected to be localized on classical regular orbits (tori), such 
as those shown in Fig. 1(b). Our numerical results indicate that {T{.\(fb) for 6 = 3 (6 = 4) 
at 7 = 73 is essentially the same as (R| for 6 = 4 (6 = 3) at 7 = 71. This "exchange" 
of Wannier functions in the passage through the degeneracy can be understood as follows. 
Sufficiently far from the degeneracy region the bands vary almost linearly as a function of 7 
(see Fig. 11 in Ref. 0) and are well approximated by a primitive semiclassical quantization 
of the two tori on which the functions {R\(pb), b = 3, 4, are localized. Near the degeneracy 



region, however, the actual band structure results from an "avoided crossing" between the 
two bands [see Fig. 12(a) in Ref. 0], leading to the exchange phenomenon. In Fig. 6(a), we 
plot the BHDs for the two bands along the symmetry line u = v for 7 = 71. These BHDs 
exhibit, essentially, the profiles of the two Wannier functions. A similar plot for 7 = 72 
(between the two degeneracies) is shown in Fig. 6(b). It is evident that the BHDs have 
changed significantly following the small variation in 7 from 71 to 72. However, Fig. 6(c) 
shows that the generalized BHD for the two bands is conserved to high accuracy under this 
variation. This BHD, associated with a GB having a total Chern index ctgb = (see Sec. 
IV), corresponds essentially to a "superposition" of the two Wannier functions. Between 
the two degeneracies, the GB can be still described by Wannier functions, associated with 
a new basis (|18D with ct„ = 0. 

As a second example, we consider a GB with (Tgb 7^ in the same model (p = 1/11). 
For 7 = 7q ^ 0.3387, a degeneracy between the bands 6 = 5 and 6 = 6 takes place. Due to 
the special symmetries of the KH model a degeneracy between bands 6 = 6 and 6 = 7 
must occur at precisely the same value of 70 but at a different value of wq. As 7 is varied 
through 7o, the Chern indices of bands 6 = 5, 7 both change from o";, = to (Xf, = —1 while 
(Jq changes from erg = 1 to ae = 3. The BHDs for all these bands are shown in Fig. 7(a) 
(for 7 < 7o) and in Fig. 7(b) (for 7 > 70). Since the bands 6 = 5, 6, 7 degenerate at the 
same value of 70, it is clear from our previous analysis that only the generalized BHD of all 
these three bands can be expected to be approximately conserved in the passage through 
the degeneracy point. In fact. Fig. 7(c) shows that this BHD, associated with a GB having 
cgb = 1, is conserved with significantly better accuracy than in the case of Fig. 6(c). 

VI. BHDs UNDER SMALL VARIATIONS OF p 

In this section, we show how to define generalized BHDs for g 7^ 1 that are continu- 
ous under small variations of p on the rationals. The topological characterization of the 
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classical-quantum correspondence in the torus is then generalized to g 7^ 1 . We shall use the 
renormalizat ion-group approach of Wilkinson [^3|,^, which was applied to the investigation 
of the spectrum of a general class of time-independent Hamiltonians on the torus. In what 
follows, we assume that the generator G of the evolution operator U = exp(— -iG) belongs 
to this class of Hamiltonians. 



We first briefly summarize the main results of Refs. |2^j2J]. Let p' = q' /p' be a rational 
number sufficiently close to p = q/p and such that p' ^ p. The p' bands for p' = q' /p' can 
then be grouped into p "clusters" of adjacent bands, where each cluster is associated in a 
natural way with a band b for p = q/p. Namely, the energy or quasienergy interval covered 
by the bands in the cluster is relatively close to that covered by band b; in addition, the 
total Chern index of the cluster is equal to cxb, 

d{b)+Nt-l 

a{C,)^ J2 (r',, = ah, (23) 

b'=d{b) 

where Cf, denotes the cluster corresponding to band b, d{b) is the label of the lowest band in 

the cluster, Nh is the number of bands in the cluster, and a^, is a Chern index for p' = q'/p'. 

The spectrum and eigenstates in cluster Cb can be approximately calculated from an effective 

Hamiltonian, Hes, obtained by properly quantizing the band function uJb{w). The effective 

scaled h in this quantization turns out to be 

Qes pq'-p'q . . 

Peff = = - — , (24 

Pes PCTb + qPb 

where pb is the integer uniquely determined from the Diophantine equation pab + qpb = 1- 
It is easy to show that the numerator and denominator in the last fraction in ( P^ are 
relatively prime integers. Eq. (|2^ ) then implies a simple formula for the number of bands 
in the cluster, Nb = Pcs, 

Nb = pob + gVft • (25) 

Now, the existence of an approximate effective Hamiltonian iYefi for a cluster means 
that the space of states in band b approximately coincides with the space of states in the 
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Nfy bands of cluster Cb- Provided p' is sufficiently close to p, any state in the cluster is 
well approximated by a linear combination of states that belong only to band b. This fact is 
expressed concisely by the statement that the projection operator for band b is approximately 
equal to that for cluster Cb- 

I I- d{b)+Nb-l ^ 

d^\^>b,^){^>b,^\^U rf^/ dW\%^^,){%^^,\, (26) 

where all the primed quantities refer to p' and A/" is a normalization constant that remains 
to be determined. Let us assume, for simplicity, that the phase-space variables R' = {u' , v') 



for p' ([«', v'\ = 2'nip') are related to the variables R = (-u, v) for p by R' = y p/ p' R. Using 
then the expressions in Eqs. and for both p and p', we immediately obtain from Eq. 
(H) that 



I . I d{b)+Nb~l ^ 



where the constant M in Eq. ( pG]) has been determined as M = l/Nf, from the normalization 
condition ([I^). Eq. ( p7D shows that the BHD Pfe(R) for band b is approximately equal to 



the generalized BHD P^^(R') for the cluster Cb- In the limit p' —>■ p, the space of the cluster 
becomes identical to that of band b [the approximate equality in Eq. is replaced by 
an equality], and P(^^(R') P;,(R). This expresses the continuity of the generalized BHD, 
associated with clusters Cb, under small variations of p on the rationals. 

The most important reference values of p are, of course, those with q = 1 (i.e., p = 
for which the topological characterization of the classical-quantum correspondence on the 



torus is well established |§-|T0[- Using the analysis in the previous paragraph, this charac- 
terization can be easily extended to rational values p' = q' /p' {q' ^ 1) sufficiently close to 
p = 1/p- Consider ffist the case of cxb = 0, with eigenstates localized on a regular classical 
orbit and exhibiting weak sensitivity to variations in the BCs. Here, the Diophantine equa- 
tion pab + qpb = 1 (Q' = 1) implies that pb = 1. Then, from Eq. (P^D, the corresponding 
cluster Cb consists simply of Nb = q' bands with total Chern index cr(Cfe) = [see Eq. (p3|)]. 
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As shown in Sec. IV, one can find for sucli a cluster (or generalized band) a new basis 
(p!8|) whose effective Chern indices all vanish. The cluster can then be viewed as exhibiting 
"weak" sensitivity to variations in the BCs, in analogy to the single band b. In addition, 
because of Eq. (0), the generalized BHD for the cluster is approximately equal to the 
BHD for band b, and is therefore localized on the same classical regular orbit. In the case 
of cTfc 7^ 0, the number A*";, of bands in the cluster is generally different from q' and varies 
with b. However, the Chern index o-{Cb) of the cluster is still equal to [Eq. (p3|)], and, 
like the single band b, the cluster exhibits "strong" sensitivity to variations in the BCs, 



in the sense that a new basis (|T8D whose effective Chern indices all vanish does not exist 
(see Sec. IV). If the BHD for band b is spread over the classically chaotic region, the same 
will be true for the generalized BHD of the cluster. The topological characterization of the 
classical-quantum correspondence on the torus is thus generalized to p' sufficiently close to 
p = 1/phj replacing single bands b with the corresponding clusters Cb- 



As an example, we consider the case of p = 1/5 for the KH model with 7 = 0.3. While 
for this value of 7 the chaotic region occupies a significant portion of the phase space (see 
Fig. 1), the Chern indices turn out to be the same as in the integrable limit (Harper model), 
i.e., o"b = for 6 = 1, 2, 4, 5 ("regular-motion" bands) and (T3 = 1 ( "chaotic- motion" band). 
Similarly, for three rational approximants p' of p, p' = 2/11, 4/21, 6/31, the Chern indices 
are the same as those in the Harper model Jl^ (however, the quasienergy spectra and states 
are, of course, completely different from those in the integrable limit). Using then formula 
(p5|), we find, for these values of p', that A'f, = q' for 6 = 1, 2, 4, 5 while = q' + 1 for 
6 = 3. Fig. 8(a) shows the BHD Ptijl) for 6 = 1 and the corresponding generalized BHDs 
Pq^(R!), calculated using Eq. (p7D, along the symmetry line u = v. Similar results for 6 = 3 
are shown in Fig. 8(b). The convergence of P^^(R') to Pf,(R) as p' approaches p is clear in 
both cases. 
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VII. CONCLUSIONS 

The framework of quantum dynamics on a toral phase space offers the possibihty of clas- 
sifying the eigenstates as "regular" or "chaotic" according to their sensitivity to variations 
in the boundary conditions (BCs), i.e., variations of the quasicoordinate w in the Brillouin 
zone (BZ). A measure of this sensitivity for a given band b is the Chern index cr^,, which 
determines the quasiperiodicity conditions (H) and satisfied by the band eigenstates. 
Provided p (the scaled h) assumes a rational value of the form 1/p {p being a sufficiently 
large integer), eigenstates well localized on regular classical orbits belong to bands with 
(Xh = 0. On the other hand, eigenstates concentrated in chaotic regions of the phase space 
usually belong to bands with cxb 7^ 0. 

By definition, at characterizes the entire band b, rather than individual eigenstates. 
The Chern-index classification into "regular" or "chaotic" refers then to bands and not to 
eigenstates. It is also well known that cr^ can be expressed as a uniform average over 

w in the BZ: 

d%^^{R) gvi>,,^(R) _ dn^^jR) gvl/,,^(R) l 

dwi dw2 dw2 dwi 

where \E'b_w(R') is the coherent-state representation of the eigenstate l^&t.w)- These observa- 
tions motivate one to introduce other global characterizations of the bands which take into 
account all the BCs only in an average sense, thus resolving the ambiguity in the choice of a 
particular BC. These additional characterizations should provide a more detailed description 
of the classical-quantum correspondence on the torus without referring to the individual, 
nonclassical BCs. 

In this paper, we have introduced and studied one such characterization, the band Husimi 
distribution (BHD), given by the uniform average of |\E'b_w(R')P over w. The BHD may be 
viewed as the representative Husimi distribution for "level" b and is expected to be closer to 
a classical probability density than |\E'b,w(R')P- This expectation was shown to be fulfilled 
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(Jh = — / dw / (i^R 

27T Jbz JtA 



in several aspects in Sec. III. The BHD concept plays an important role in extending the 
Chern-index characterization of the classical-quantum correspondence to rational values of 
p with numerators larger than 1 (see Sec. VI). In this case, a generalized version of the 
BHD has to be introduced (see Sec. IV), given by the average of the BHDs associated with 
a cluster of adjacent bands for p' = q'/p' ^ 1/p. 



Smoothing a quantum probability distribution over a range of energy levels is impor- 
tant in the theoretical study of scars using the semiclassical periodic-orbit theory In 
the case of the BHD, this smoothing is performed over the continuous range of one band, 
corresponding essentially to a single level in the framework of a toral phase space. Using 



the adaptation of periodic-orbit theory to this framework ||T2[, it may be possible to achieve 
a better understanding of the nature of BHDs and generalized BHDs in the semiclassical 
limit. 



An important question is to what extent information about the individual band eigen- 
states can be recovered from the BHD. We recall here the well known fact in solid-state 
physics that all the Bloch eigenstates in a band can be simply reproduced from the corre- 
sponding Wannier function. Thus, in this case, no information about the individual eigen- 
states is lost by averaging over them (the average being the Wannier function). Rel. (|T5|), 
valid only for o";, = 0, is analogous to the formula expressing the Bloch eigenstates in terms 



of the Wannier function. In this case, Rel. (|T^) shows that the BHD is closely related to 
the absolute value squared of a Wannier function. This relation may be extended to o";, 7^ 
using results from Ref. |Q. It was also shown recently |Q that the Wigner functions of 
all of the band eigenstates can be reproduced from the Wigner analog of the BHD. On the 
basis of these observations, we expect that it should be possible to recover, at least partially, 
relevant information about the individual eigenstates from the BHDs. This will be investi- 
gated in future works. 
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FIGURES 



FIG. 1. Classical Poincare maps of typical orbits of the kicked-Harper (KH) model (|1]) 
for different values of the nonintegrability parameter 7 = Ar/(27r/): (a) 7 = 0.001 (nearly- 
integrable regime), (b) 7 = 0.26 (mixed regime). Notice, for future reference, the island 
chains of periods 6 and 8 surrounding the central elliptic point, (c) 7 = 0.56. The chaotic 
region occupies a large fraction of the phase space, (d) 7 = 0.95 (strongly-chaotic regime). 

FIG. 2. Density plots of Husimi distributions |\I'b,w(R')P for the KH model with pa- 
rameters 7 = 0.26 [compare with Fig. 1(b)] and p = 1/11. The central band b = 6 (with 
Chern index erg = 1) is considered for different values of w in the BZ: (a) w = (0, 0), (b) 
w = (vr/ll, 0), (c) w = (0, vr/11), and (d) w = (tt/II, vr/ll). In these plots, as well 
as in Figs. 3-5, we use a power-law density scale (with power 1/3) and ten gray tones, 
with darker tones corresponding to higher values of the distribution. In cases (b) and (c), 
|^'b^w(R')P is concentrated on one of the two hyperbolic fixed points, while in cases (a) and 
(d) it is concentrated on both points. 

FIG. 3. Similar to Fig. 2, but for band 6 = 5 with Chern index = 0. Despite the fact 
that (J5 = 0, we observe strong sensitivity to variations in w, as in the case of Fig. 2. In plot 
(a) [w = (0, 0)], |\E';,,w(R)P is concentrated on the island chain of period 8 [see Fig. 1(b)]. 
Notice that the localization regions in the case of plots (b) and (c) are reversed relative to 
the corresponding plots in Fig. 2. Plot (d), on the other hand, is qualitatively similar to 
Fig. 2(d). 

FIG. 4. Band Husimi distribution (BHD) for the case considered in Fig. 2 (6 = 6). The 
BHD was calculated by averaging 1^5^(1^)1^ over 20 x 20 values of w, uniformly distributed 
in the BZ. This BHD, with a localization region similar to that of Figs. 2(a) and 2(d), is 
close to a classical probability distribution for the stochastic layer in Fig. 1(b) (see text). 
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FIG. 5. Similar to Fig. 4, but for the case considered in Fig. 3 {b — 5). The BHD appears 
to be concentrated on four of the six islands of an island chain of period 6, surrounding the 
main island chain of this period in Fig. 1(b). The localization region for this BHD is 
representative for band b — 5 but is qualitatively different from that of |^b,w(R-)P for all 
the special values of w considered in Fig. 3. 

FIG. 6. (a) BHDs along the symmetry hne u — v ior bands b — 3 (solid hne) and 
b = 4 (dashed line) in the case of the KH model with parameters 7 = 0.26 and p = 1/11. 
The BHDs were calculated by averaging |^'b^w(R-)P over 100 x 100 values of w uniformly 
distributed in the BZ. (b) Similar to (a), but for 7 — 0.2645 > 70 0.264, where 70 is 
a degeneracy point for bands 6 = 3, 4. Notice the dramatic change of the BHDs relative 
to those of case (a), despite the small variation in 7 through the degeneracy point, (c) 
Generalized BHD for bands b — 3, 4, obtained by averaging the results for the two bands 
in case (a) (solid line) and in case (b) (dashed line). 

FIG. 7. (a) Similar to Fig. 6(a), but for the bands 6 = 5 (solid line), 6 = 6 (dashed line), 
and b — 7 (dotted line) in the case of 7 = 0.3385 < 70 f=:i 0.3387, where 70 is a degeneracy 
point of all the three bands, (b) Similar to (a), but for 7 = 0.339 > 70. (c) Generahzed 
BHD for bands 6 = 5, 6, 7, obtained by averaging the results for the three bands in case (a) 
(solid line) and in case (b) (dashed line, essentially indistinguishable from the solid line). 

FIG. 8. BHDs along the symmetry line u — v ior bands and clusters of bands in the KH 
model with nonintegrability parameter 7 = 0.3: (a) The upper curve (see a magnification in 
the inset) is the BHD for band 6 = 1 in the case of p = 1/5. The other three curves are the 
generahzed BHDs for the corresponding cluster C5 in the cases ofp = 6/31, 4/21, 2/11, in 
order of descending curves (see also inset), (b) Similar to (a), but for the central band 6 = 3. 
In both cases, the BHDs were calculated by averaging |^'b^w(R')P over 400 x 400 values of 
w uniformly distributed in the BZ. 
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